Ising transition in the two-dimensional quantum Ji — J2 Heisenberg model 
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We study the thermodynamics of the spin-S two-dimensional quantum Heisenberg antiferromag- 
net on the square lattice with nearest (Ji) and next-nearest (J2) neighbor couplings in its coUinear 
phase (J2/J1 > 0.5), using the pure-quantum self-consistent harmonic approximation. Our results 
show the persistence of a finite-temperature Ising phase transition for every value of the spin, pro- 
vided that the ratio J2/J1 is greater than a critical value corresponding to the onset of coUinear 
long-range order at zero temperature. We also calculate the spin- and temperature-dependence of 
the coUinear susceptibility and correlation length, and we discuss our results in light of the experi- 
ments on Li2VOSi04 and related compounds. 
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The study of frustrated quantum spin systems is one 
of the most challenging and exciting topics in theoret- 
ical magnetism. A very extensively investigated, yet 
largely debated model is the so-called J1 — J2 Heisen- 
berg model with competing antiferromagnetic couplings 
(Ji, J2 > 0) between nearest-neighbors (n.n) and next- 
nearest-neighbors (n.n.n.) 



H = Ji ^ Si ■ Sj + J2 ^ Si ■ S. 
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where Si are spin-5 operators on a periodic lattice with 
N ^ L X L sites; hereafter a = J2/J1 defines the frustra- 
tion ratio. 

In the classical limit (S 00), the minimum energy 
configuration has conventional Neel order with magnetic 
wave vector Q = (tt, tt) for a < 0.5. Instead, for a > 0.5, 
the antiferromagnetic order is established independently 
on the two sublattices, with the two (staggered) magneti- 
zations free to rotate with respect to each other. Among 
this degenerate manifold, two families of coUinear states, 
with pitch vectors Q = (7r,0) or (0,7r), are selected by 
an order-by-disorder mechanism as soon as thermal or 
quantum fluctuations are taken into account. As a re- 
sult, for a > 0.5 the classical ground state breaks not 
only the spin rotational and translational invariance of 
the Hamiltonian - as the conventional Neel phase - but 
also its invariance under tt/2 lattice rotations, the result- 
ing degeneracy corresponding to the group 0(3) x Z2. 
Remarkably, the additional discrete Z2 symmetry can in 
principle be broken at finite temperatures without violat- 
ing the Mermin- Wagner theorem , which applies in two 
dimensions only to continuous ones. On this basis, in a 
seminal paper Chandra, Coleman and Larkin (CCL) 
proposed that the two-dimensional Ji — J2 model could 
sustain a Ising phase transition at finite-temperature, 
with an order parameter directly related to the Z2 de- 
gree of freedom induced by frustration. They also pro- 



vided quantitative estimates of the critical temperatures 
in the large-a limit for both the classical and the quan- 
tum cases. 

While the CCL transition in the classical model has 
been recently established by an extensive Monte Carlo 
(MC) study 01, the occurrence of a low-temperature 
phase with a discrete broken symmetry in the quantum 
case is still a subject of debate. Besides its own theoreti- 
cal interest, this issue has become particularly important 
in connection with the discovery of three vanadate com- 
pounds (Li2VOSi04, Li2VOGe04, and VOM0O4) whose 
relevant magnetic interactions involve nearest and next- 
nearest spin-1/2 ions on weakly coupled stacked 
planes 4,:5|. In particular, NMR and /iSR measurements 
on Li2VOSi04 4'| indicate the occurrence of a transi- 
tion to a low-temperature phase with coUinear order at 
~ 2.8 K. The selection of the coUinear ground state 
suggests the CCL mechanism as possible underlying ex- 
planation. Unfortunately, a clear experimental and theo- 
retical picture is still elusive: in the experiments with 
vanadate compounds, structural distortions, interlayer 
and anisotropy effects are likely to come into play 4], 
and on the other hand the theoretical investigation can- 
not rely on the insight provided by quantum Monte Carlo 
methods as their reliability in presence of frustration is 
strongly limited by the infamous sign problem 0. 

The necessary condition for the CCL transition to take 
place is the presence of coUinear order at zero tempera- 
ture. To this respect, the existent theoretical results 
point towards a collinear-ground state for frustration ra- 
tios a > ttc, with the critical value, ac, increasing as the 
value of the spin decreases (for S — 1/2, — 0.6). Be- 
low this value quantum fluctuations seem to be strong 
enough to stabilize a low-temperature phase with short- 
range magnetic correlations |7| , but above it a CCL tran- 
sition is in principle possible. However, this possibility 
has been recently challenged by various high-temperature 
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expansion studies H, 0, which were not able to de- 
tect any evidence of a finite-temperature transition for 
S — 1/2. On this basis, Singh et al. have recently 
proposed a scenario where, due to quantum fluctuations, 
the broken lattice symmetry of the coUinear ground state 
would be restored at any non-zero temperature. 

In this paper, we present a complete study of the ther- 
modynamic properties of the the quantum J1 — J2 model 
in its collinear phase, obtained within an effective Hamil- 
tonian approach: the pure-quantum self- consistent har- 
monic approximation (PQSCHA) This approach, 
based on to the path-integral formalism, allows one to 
separate the classical from the pure-quantum contribu- 
tion to the thermodynamics of the system. Both the 
classical physics and the purely quantum linear effects are 
exactly described within the PQSCHA at any tempera- 
ture, while the purely quantum non-linear contributions 
are treated within a self-consistent harmonic approxima- 
tion. This feature makes the PQSCHA a valid tool to 
investigate the effects of quantum fluctuations on phase 
transitions ~ like the CCL one - whose character is essen- 
tially classical. In particular, this approach has been suc- 
cessfully applied to a variety of spin systems displaying 
Kosterlitz-Thouless and/or Ising critical behaviors, pro- 
viding reliable estimates of the transition temperatures 
even for = 1/2 0- 

Within the PQSCHA framework, the thermodynamics 
of a quantum system is rephrased in terms of a classical 
effective Hamiltonian with renormalized parameters de- 
pending on the spin value, temperature, and frustration. 
The derivation of the effective Hamiltonian for the J1 — J2 
model closely follows the steps shown in Ref. (Hj. The 
only detail which is worth mentioning here is that, in 
this case, the calculation of quantum renormalizations, 
involving a harmonic expansion around one of the two 
families of collinear states, gives rise in general to so- 
lutions with an explicitly broken symmetry under tt/2 
lattice rotations. However, it is possible to show that, 
to 0{1/S), the effective Hamiltonian can be recast in a 
form preserving all the symmetries of the original model, 
and that reads (except for uniform terms): 
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where Si are classical vectors of length 1, S = S ■ 
the effective spin length "ll], and jf = {91 -\- el)ejJi/2, 
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— O2J2, are the quantum-renormalized exchange 



integrals, with spin-, temperature- and frustration- 
dependent renormalization parameters 6x, Oy and ^2- 
These, e.g. referred to the ground state with Q = (tt, 0), 
are given by 0^ = 1 — 1^01/2 {a = x, y, 2) where the coef- 
ficients 
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FIG. 1: Renormalized critical temperature of the CCL tran- 
sition for various values of the spin evaluated using Eq.i 
(non-solid lines). Classical data (•) are taken from Ref. 
The solid lines on the right are the CCL and the CCL* pre- 
diction for the classical case (see text). The arrow marks the 
boundary of the non- magnetic phase for S — 1/2 p^ . 



are self-consistently determined with 
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COS kx cos ky, T is the temperature, and the renormalized 
dispersion relation is ujf, = '^k'^k- Interestingly, at T = 
the PQSCHA turns out to be equivalent to the mod- 
ified spin-wave theory 0, which is expected to be a 
faithful representation of the low-energy properties in the 
collinear phase. In particular, the spin-wave dispersion 
relation ujk has been recently shown to be in remarkable 
quantitative agreement with series expansion results in 
the entire range of momenta for a > ac, and it is there- 
fore expected to provide an accurate description of the 
low-energy excitations of the model . 

The occurrence of the CCL transition in the quantum 
case can be directly addressed within our approach by 
calculating the critical temperatures as functions of the 
spin and of the frustration ratio: a zero value of the crit- 
ical temperature, or the breakdown of the self-consistent 
harmonic treatment of quantum fluctuations, would sig- 
nal a possible absence of the phase transition. In par- 
ticular, using a simple scaling argument the critical tem- 
peratures in the quantum case Tc{S, a) can be related to 
those of the classical model Tc'^^\a) through the follow- 
ing self-consistent relation |3| 

n{S, a) = jfin, S, a) Ti^')(a^ff (Te, S, a)) , (5) 
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where jf = jfS^/Ji and a"^ = Jf/Jf- The 
classical transition temperature, Tjf'^\a) is accurately 
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FIG. 3: PQSCHA results for the specific heat in the cases 
FIG. 2: PQSCHA resuhs for the specific heat (upper panel) a = 1, S = 1, 1/2. In the inset: scaling of the specific heat at 
and the order parameter a (lower panel) as function of tern- ^he peak temperature ksT/JiS^ = 0.36 for S = 1. 
perature for various spins, in the case a = 0.65. 



known through extensive MC simulations for a < 2; 
it vanishes for a — + 1/2 and grows more or less lin- 
early for a > 1. Beyond a ~ 2 the determination 
of the critical temperatures becomes troublesome due 
to severe finite-size effects related to the width of the 
domain walls between domains with Q = (tt, 0) and 
Q = (0,7r) 0. However, for large values of the frus- 
tration ratio the classical CCL estimate of the critical 
temperature, = 0.768J2/[1 0.135 ln(J2/Ji)], is ex- 
pected to be reliable 0, 0] . This vanishes logarithmically 
in the limit J2/J1 — * 00, corresponding to two decou- 
pled unfrustrated Heisenberg systems. In order to rep- 
resent the whole interval of a G [1/2, 00) in Fig. ^ we 
have plotted both the MC and the CCL estimates of the 
classical critical temperatures as a function of a/(l + a). 
The mismatch between the MC and CCL predictions is 
a minor flaw that can be easily accounted for, and cor- 
rected by slightly modifying CCL's criterion for the de- 
termination of the transition temperature as explained 
in Ref. This gives rise to the curve marked in Fig.Q] 
as CCL*. Using the latter estimate, and the classical 
MC results of Ref. Q , the transition temperatures in the 
quantum case can be determined by numerically solving 
Eq. inj. As shown in Fig. the critical temperatures 
decrease as the spin decreases, due to the enhancement 
of quantum fluctuations. Remarkably, while for large a 
the transition temperature vanishes for a 00 for any 
value of the spin, in the opposite limit the critical tem- 
peratures vanish approaching a critical value ac > 0.5 
that increases as S decreases, thus confirming the ex- 
istence of a non-magnetic phase in the regime of high 
frustration. In particular, for S = 1/2, ac — 0.6 in agree- 
ment with the previous estimates of the zero-temperature 
quantum critical point For large a, the PQSCHA re- 
sults turn out to be consistent with the CCL estimates 
for the quantum S' = 1/2 system |^|^. For intermediate 
values of a, the transition temperatures remain finite for 



any spin value. However for S = 1/2, though sizeable 
{Tc/J2 — 0.2), the critical temperatures are more than 
an order of magnitude smaller than those considered in 
Ref. ;9j in the discussion of the high-temperature expan- 
sions results leading to the proposal of a T = critical 
scenario. On the contrary, a finite-temperature phase 
transition at the critical temperatures we estimate is in 
fact consistent with the numerical results of Ref. |^ . 

The study of the thermodynamics of the effective clas- 
sical model ^ provides a deeper analysis of the CCL 
transition. To this end, we have performed classical MC 
simulations on the effective Hamiltonian on L x L lat- 
tices, with L up to 300. The PQSCHA expression for 
the order parameter associated to the CCL transition Q 
is CT = 6*1 {0l + el){\{si - S3) • (S2 - S4)|)eff/8, where 
(1,2,3,4) are the sites of on an elementary plaquette of 
the lattice with counterclockwise numbering, and (...)eff 
is the statistical average associated to the effective classi- 
cal Hamiltonian Another quantity bearing clear sig- 
natures of the transition is the specific heat, calculated 
here as the numeric derivative of the internal energy per 
spin, c = du/dT, where u = {nes)/N - JiS^iO^ - O^). 
These quantities are shown for a = 0.65 in Fig. |5| The 
Ising transition is clearly marked by a (finite-size) peak 
in the specific heat in correspondence to the temperature 
where the Ising order parameter vanishes. By decreasing 
the spin, the saturated values of a decrease as well, due 
to the increased quantum fluctuations, and so does the 
peak of the specific heat, making the transition in the 
quantum case generally weaker than in the classical one. 
When a is increased, the signatures of the transition are 
even more dramatically suppressed, due to the mentioned 
effects related to the width of the domain walls . The 
Heisenberg features become more and more important 
with respect to the frustration effects, in fact, as shown 
in Fig. 13 already for a = 1, the Ising peak of the specific 
heat tends to be masked by the broad maximum remi- 
niscent of the unfrustrated limit a ^ 00, and very large 
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FIG. 4: Upper panel: 5 = 1/2 staggered susceptibility, with 
Q = (tt, 0) and (tt, tt) in the frustrated (open symbols) and 
unfrustrated cases, respectively. Lower panel: 5 = 1/2 corre- 
lation len gth . The data for the unfrustrated cases are taken 
from Ref. I 111 and partly from Ref. llTl 

values of L are required in order to resolve the logarith- 
mic divergence (inset). Similarly, the absence of critical 
features in the recent numerical calculations of the spe- 
cific heat |0| can be traced back to the limited range 
of the correlators generated by the high-temperature ex- 
pansion [l5| . 

In order to investigate the effect of frustration on the 
antiferromagnetic correlations, we have also calculated 
the collinear susceptibility 

X(Q) = + ^ E ^^''■'^(^0 ■ 'rU , (6) 

with Q = (tTjO) or (0, tt), and the corresponding correla- 
tion length (through the second moment estimator ^Td\). 
The results for the case S = 1/2 for different values of 
a are shown in Fig. ^ We observe that both quantities 
are strongly suppressed by frustration for a ^ 1, but for 
larger a the results are extremely close to those in the un- 
frustrated Heisenberg model, corresponding to a ^ oo. 
Therefore neutron scattering experiments on Li2VOSi04 
and related compounds with a large a are not likely to 
show any strong frustration effects in the magnetic cor- 
relations, as also indicated by preliminary experimental 
results 

Finally, we compare the estimated critical temperature 
for S = 1/2 with the transition temperature, T„ = 2.8 K, 
observed in Li2VOSi04. For this compound, various ra- 
tios J2/J1 have been estimated, ranging from J2/J1 = 1-1 
with Ji = 3.9 K to J2/J1 = 4.76 with Ji = 1.25 K 1^. 
Using these estimates, we get Tc « 1.01 K and 1.46 K, 
respectively, which are well below the transition to three- 
dimensional collinear order observed at T., ~ 2.8 K. The 



CCL critical behavior is therefore not detectable in this 
compound 0. In general, due to the weak nature of the 
transition for J2/J1 ^ 1, previously discussed, the obser- 
vation of the CCL transition would require realizations 
of a Ji — J2 model with a smaller value of the ratio J2/ Ji- 

In conclusion, using the pure-quantum self-consistent 
harmonic approximation, we have provided a complete 
and consistent picture of the collinear phase of the 2D 
quantum Ji — J2 antiferromagnet. Our results indicate 
that the finite temperature transition predicted by Chan- 
dra, Coleman, and Larkin i2| persists down to 5* ~ 1/2 
provided that the ratio J2/J1 is greater than a critical 
value corresponding to the stabilization of collinear long- 
range order in the ground state. We believe that our find- 
ings will be a valid reference point for future experimental 
investigations on Li2VOSi04 and related compounds. 
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